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This paper presents a non-tubular heat exchanger for use in a Stirling solar engine micro co-generation 
unit that is being developed by the University of Oviedo and the technological research centre Tekniker 
Foundation. The engine has been designed using scaling criteria, which imply a case of relaxed dynamic 
similarity, that is to say, without preserving the equality among all of the prototype's and the derivative’s 
dimensionless groups. The geometrical characteristics of the heater are described and the variables that 
can influence the pressure drop and heat transfer in the heater are identified. The correlations for the 
friction factor and Stanton number have been measured under steady flow conditions to analyse the 
feasibility of the experimental heater. The Reynolds analogy is not fulfilled, and the compressibility 
effects are negligible. The comparison between the correlation predictions for the non-tubular heater 
and slotted or tubular heaters of well-known prototypes shows possible applications for heaters that 
reach high Reynolds numbers and that encounter practical problems associated with the use of bundles 
of tubes. 

© 2012 Elsevier Ltd. All rights reserved. 


1. Introduction 

Several countries are currently investigating distributed power 
generation systems fed by renewable energy to alleviate the satu¬ 
ration of supply networks in developed areas, as well as to provide 
an alternative for the electrification of rural districts in developing 
countries [lj. 

Dish-Stirling systems are currently a commercial reality [2], 
with the highest efficiency in solar conversion at power levels that 
are interesting for distributed power generation [3]. 

The engines in these systems have been adapted from proto¬ 
types originally designed for other applications. The heaters in 
those engines are variations of tubular heater heads located at the 
focal plane of a paraboloid solar concentrator. 

The University of Oviedo and the technological research centre 
Tekniker Foundation, from Eibar, Spain, are engaged in a joint venture 
to develop a Stirling solar engine micro co-generation unit. The 
engine is currently being tested and is expected to serve low power 
applications, for which Stirling systems can be competitive [4]. 

The prototype has been designed to use two possible types of 
heaters, one of which has non-tubular geometry. The present paper 
focuses on the geometrical and physical features of a non-tubular 
Stirling engine heater, which could be interesting from both the 
design and performance viewpoints. 
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2. Characteristics of usual Stirling engine heaters 

Philips’ research has remarkably influenced the evolution of 
Stirling engine heater design. Since the fifties, the usual approach 
has consisted in slotted or tubular heaters supplied by combustion 
heating systems. Those tubes encountered problems related to 
their design, manufacturing and maintenance, such as the lack of 
thermal distribution uniformity, the risks of corrosion and over¬ 
heating, the cost of special high temperature and high pressure 
resistant alloys, the combustion products fouling, etc. The tubes 
were usually joined to the cylinders by brazing either in hydrogen 
or vacuum atmosphere because of their inaccessibility for welding. 
A typical example is the helical heater-tube cage of the 4-215 
engine tested in the Ford Torino around 1975. The tubes were made 
of Multimet with A1SI-310 fins arranged on the outside branches. 
Around 1976, Meijer asserted that the complexity of the heater 
tubes made of exotic materials must be reduced to make them 
price-competitive [5], 

However, substituting or modifying a typical tubular heater 
must preserve the flow symmetry as much as possible without 
penalising the dead volume or other parameters that influence the 
working fluid friction and heat transfer. The case of the SOLO V160 
alpha-type co-generation engine is an example of an adaptation to 
solar applications. Since the nineties, SOLO has developed 
improved units with various heaters consisting of tubes brazed to 
tower manifolds. These heaters have accumulated thousands of 
operation hours with direct solar reception in Dish-Stirling units 
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Nomenclature 

A wxe wetted area of heater, m 2 

A xxr , cross-sectional area of heater, m 2 

c p specific heat of gas at constant pressure, J/(kg K) 

c v specific heat of gas at constant volume, J/(l<g 1<) 

c w specific heat of wall material, J/(kg K) 

Cf instantaneous friction factor =N eu (RhelL e ) 

h convective heat transfer coefficient, W/(m 2 K) 

k thermal conductivity of the gas, W/(m K) 

L e characteristic heater length, m 

m mass flow rate, kg/s 

n s engine frequency, rev/s 

N eu instantaneous Euler number 

= a P/Qp^ 2 ) = 2 ( A P/P) N ma 

Nf 0 instantaneous, spatially-averag ed Fo urier 

number =a w /uR he = a w /Rh e \JR k re f ■ 1 /N ma 
N ma instantaneous, spatially-averaged, isothermal Mach 

number =u/RT ref = m^RT^f/pA^e 

N ma characteristic Mach number = n s Vsw 3 / y/RTc 
N p characteristic pressure number =N S gNma 

N pr Prantdl number =pc p H< 

N re instantaneous, spatially-averaged Reynolds 

number =4 puR he /p = 4mR he /pA xxe 
N sg instantaneous, spatially-averaged Stirling 

number =NmalN re 

Nsc characteristic Stirling number =p m lpn s 

N st instantaneous, spatially-averaged Stanton 

number =h/puc p = hA xxe /mc p 
N T cr instantaneous thermal capacity 

ratio =p w c w /pc p = p w c w T ref /p-(y - l)/y 
p pressure, Pa 

p m mean pressure, Pa 

P B brake power, W 

Pind indicated power, W 

R specific gas constant, J/(kg K) 

r he local hydraulic radius of heater, m 


Rhe characteristic hydraulic radius of heater, m =4Vd Xe IA wxe 

t time, s 

T, T g gas temperature, K 

T w wall temperature, K 

u gas velocity, m/s 

u mean gas velocity, m/s 

V ie dead volume of expansion space, m 3 

Vdxe dead volume of heater, m 3 

V cylinder volume, m 3 

V sw swept volume, m 3 

W 0 quasi-static work per cycle, J 

x longitudinal coordinate along the heater length, m 

cixxe cross-sectional area ratio =A xxe lvM? 

a w thermal diffusivity of wall material, m 2 /s 

a wxe wetted area ratio =A wxe /V 2 ^ 

/? drag coefficient in Forchheimer’s law, m _1 

Ai,...,A fc dimensionless geometrical parameters, including 

those characteristic of the drive mechanism 
4> coefficient of indicated power losses 

y adiabatic coefficient 

k swept volume ratio =V C IV E 

k permeability coefficient in Forchheimer’s law, m 2 

Ahe hydraulic radius ratio =r/, e /V'i^ 

p viscosity, Pa s 

Pde dimensionless dead volume of expansion space =Vdel 

V sw 

Pdxe dimensionless dead volume of heater =Vd Xe IV sw 

W coefficient of indicated power losses 

p density of gas, kg/m 3 

p w density of wall material, kg/m 3 

t temperature ratio =T C /T £ 

Cind dimensionless indicated power =Pmdl(PmV sw ns ) 

Co quasi-static dimensionless work per cycle =W 0 / 

C PmVsw ) 

Subscripts 

c, C compression space 

e, £ expansion space 


[6], However, tower manifolds penalise dead volume and flow 
symmetry, which suggests that other approaches might increase 
both the reliability and the high efficiency of the SOLO units. 

The flow symmetry is more adequate in the heater cage designed 
by Mahkamov et al. [7] for a 3 kW Stirling engine, which combines 
natural gas and solar energy in remote sites of Central Asia. The 
combustion chamber provides a cavity that surrounds the heater 
and has an input window at its top for the concentrated solar irra¬ 
diation. The heater cage consists of two collectors and two sets of 
semi-circular rings made from stainless steel tubes. A different 
approach uses a direct solar heater-tube receiver placed at the 
bottom of a cylindrical cavity in the multi-cylinder United Stirling 4- 
95 engine, which holds the solar-electric efficiency world record [8]. 

To reduce the geometrical limitations when adapting heater- 
tube cages to solar conversion, heat-pipe receivers have been 
suggested as intermediate heat transfer systems. The STM4-120 
prototype adapted by Stirling Thermal Motors is an example 
technology that allocates a solar dish heat-pipe receiver in multi¬ 
cylinder engines at high power ranges. At a low power range, the 
NALSEM 500 unit, originally conceived to provide solar power for 
space stations, uses an interface of sodium heat pipes to uniformly 
transfer heat from a cylindrical solar cavity receiver to the heater 


head [9]. However, intermediate heat transfer systems have 
disadvantages for both the system’s efficiency and costs, particu¬ 
larly for low power, terrestrial applications. 

Additionally, despite the fact that the most efficient solar Stirling 
engines have been developed by Philips licensed companies, the 
design methods, comprehensive specifications and test results 
related to their prototypes have not been published. Independent 
authors, like Tew et al. [10], Tsuchiya et al. [11], Rice and Hennes 
[12], Lee et al. [13], and Organ et al. [14], have studied the gas 
dynamics and heat transfer of tubular or slotted Stirling heaters, but 
there are no current friction factors or Stanton or Nusselt number 
correlations that are corroborated enough for not-fully-developed, 
bidirectional flow with a variable mass rate. 

The publications related to other geometries may be considered 
exceptional and incomplete. For example, Isshiki et al. [15] showed 
experimental results for an engine with a tubeless heater, but the 
heater performance characteristics were not reported. 

In short, the Stirling engine heaters designed for solar applica¬ 
tions usually have tubular arrangements whose development level 
does not seem definitive; meanwhile, innovations are hindered by 
the lack of knowledge about the performance of new geometrical 
concepts. 
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3. Geometrical description of the experimental heater 

Similarity criteria previously published have been applied to the 
preliminary design of the Stirling engine developed by the 
University of Oviedo and the Tekniker Foundation [16—21]. The 
engine was scaled from the Philips M102C beta prototype [22], 
using air as the working fluid and a ‘bell’-crank mechanism. The 
mean pressure and engine speed were assumed to have more 
influence on the engine’s cost than the swept volume. Thus, 
a length scale factor of 1.6 was assumed, which should lead to 
a more powerful, slower and less pressurised engine using the 
same working fluid [23,24]. Table 1 shows the main characteristics 
and operating points for both the prototype and the derivative, 
assuming strict dynamic similarity, that is to say, the equality 
between all the dimensionless groups of the prototype and the 
derivative. 

Scaling-based design methods provide criteria to adapt an 
engine for an application different from its original one. However, 
heater innovation is difficult without using out-to-scale geomet¬ 
rical parameters. Thus, the concept of ‘relaxed dynamic similarity’, 
introduced by Organ [16,17], should be taken into account. This 
concept produces an interesting derivative despite not maintaining 
the equality between all of the dimensionless groups of the 
prototype and the derivative. The practical rules for applying this 
procedure have been published [24]. For example, a model might 
be interesting with a heater whose hydraulic radius-to-length ratio 
is higher than the prototype value, which means the lack of 
geometric similarity, provided that the approximate equality 
among the rest of dimensionless groups, leads to an acceptable 
power. 

The strict similarity was only maintained for the cylinder 
diameters, the piston strokes, and the geometrical characteristics of 
the slotted cooler, as well as for the wire diameter, porosity, free 
flow area and length of the regenerator because the compactness 
and low speed operation may justify the target of an alpha-type 
derivative with the Ross yoke mechanism. The combination of 
these design criteria led to an out-to-scale swept volume of 
V'sw = 341.82 cc, which implies a case of relaxed similarity. Thus, the 
validity of the scaling procedure can only be accepted if the lack of 
strict similarity in several components does not lead to character¬ 
istic dimensionless parameters that would be unfavourable for the 
derivative operation. 

Fig. 1 shows a schematic section of the experimental engine. The 
buffer space is enclosed at the bottom by a casing, which has been 
removed in the figure to show the drive mechanism and the elec¬ 
trical generator. The heater was manufactured using die-sinking 
electrical discharge machining (EDM) to form the required shape 
in the fiat disc of stainless steel (shown at the top level of the figure) 
with a three-dimensional electrode of graphite. 

The heater geometry consists of a flat surface intended to receive 
direct solar irradiation from a conventional optical concentrator, for 
example, a parabolic dish or a Fresnel lens-based system, that will 
then transmit it to the bases of almost one thousand 10-mm long and 
2-mm diameter rods, around which the working gas circulates. The 
rods are arranged in a staggered manner with a uniform diagonal 
pitch on an oval area perpendicular to the cylinder axes, with ends 


Table 1 

The main design parameters of the prototype and the derivative under strict simi¬ 
larity conditions. 



Te 

Tc 

V E 

Psw 

Pm 

n s 

Pind 

Pb 


(°C) 

(°C) 

(cm 3 ) 

(cm 3 ) 

(MPa) 

(rpm) 

(W) 

(W) 

M102C 

600 

60 

59.4 

64.2 

1.103 

-I860 

-375 

-163 

Derivative 

600 

60 

243.3 

262.8 

0.69 

-1160 

-600 

-260 
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Fig. 1 . A picture of the experimental engine. 


that are connected to the hot cylinder and the regenerator casing, as 
shown in Figs. 2 and 3. 

Table 2 allows the geometrical characteristics of the experi¬ 
mental heater to be compared with those of the prototype heater 
and the ones of a derivative designed under strict similarity 
conditions. The geometry of the experimental heater requires 
a particular definition for the heater length L e , which has been 
defined as the distance between longitudinal axes of the cylinders. 
Thus, L e is a specimen of ‘equivalent’ length that must be taken 
into account to accurately interpret the meaning of the cross- 
sectional area A xxe and those of the non-dimensional ratios Rj, e / 
Le and tr XX e- 

Dynamic similarity requires that the Stirling engine heater 
prototype and its derivative have the same values of a set of 
dimensionless parameters, namely, three independent values 
related to geometry (for example, Hd xe , a wxe and Aj, e ) and those 
related to the working fluid and operating conditions (t, y, N p and 
Nma)- 

The distinction between the scope of Thermodynamics and 
other physical theories related to Gas Dynamics and Heat Transfer 
[25] is important when dimensional analysis techniques are 
applied. Such a distinction refers to the selection of both the 
dimensional basis and the physical quantities that influence the 
problem. The thermodynamic similarity is influenced by the 
temperature, mean pressure and volume parameters; therefore, the 
quasi-static dimensionless work per cycle may be expressed as 
follows [19,24]: 

?0 = /(t)A],..., A^, Mde* Mdn Mdci Mdcc) 

To take into account the features related to Gas Dynamics and 
Heat Transfer, i.e„ the friction factor and Stanton or Nusselt 
numbers, the rest of the geometrical parameters, physical proper¬ 
ties of both the working fluid and the regenerator material, and 
operating engine speed, must be added to the list of influencing 
quantities, which allows the following expression to be obtained 
for the dimensionless indicated power [19,24]: 

Cind = /(r, K- A]..... A/ ( . tk/ce■ Ude ■ Mdr■ Mdc■ Mdcc> &xe- Ot-xr, ®xc, /-hr- 
^ he i y^N a ,N TC R,Np,N MA ) 

Prieto et al. [20] have proposed the following approach that 
explicitly states the influence of the engine speed on the dimen¬ 
sionless indicated power of kinematic Stirling engines: 
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Fig. 2. A conceptual sketch of the new heater. 


?ind = ?0 - @ N MA ~ VN 2 ma (1) 

Thus, in this equation, £o is a thermodynamic concept that 
represents the theoretical prediction, assuming an ideal cycle 
without any losses caused by leakage, thermal or mechanical irre¬ 
versibilities, while the factors $ and ¥ of indicated power losses 
characterise irreversibilities inherent to working gas friction and 
heat transfer. 

The new heater’s dead volume is more than twice the value of 
the strictly similar heater based on thermodynamic similarity. This 
disadvantage has been compensated for by reducing the dead 
volume in both transition spaces adjacent to the expansion cylinder 
and the regenerator. Therefore, the value ~ 0.302 was computed 
by quasi-static simulation for the relaxed derivative, which is 
slightly higher than the value £o = 0.294 obtained for the proto¬ 
type, with a different engine configuration and drive mechanism. 

Table 2 shows that the experimental heater ratio RhelL e is more 
favourable with regards to the pressure losses caused by the 
working gas friction, at 41% higher than the strict derivative value. 
Conversely, the heat exchanging capability of the experimental 
heater might be compromised because its wetted area ratio ct wxe is 
19% lower than the strict derivative’s value. This handicap is caused 
by the high swept volume derived from the new engine configu¬ 
ration and drive mechanism, despite the fact that the wetted 
surface is almost identical for both cases. 



Fig. 3. The rods of the experimental heater at the ends of the hot cylinder and the 
regenerator casing. 


Thus, the geometrical design of the engine seems adequate 
based on the quasi-static performance; however, before accepting 
the feasibility of the new heater, additional analyses are required to 
evaluate the influence of non-geometrical variables on indicated 
power losses. 

4. Friction and heat transfer characteristics of the 
experimental heater 


Assuming a one-dimensional model, the local instantaneous 
pressure drop through an elemental length of the experimental 
heater can be expressed as a function of: the pressure itself, the 
geometrical parameters, the fluid physical properties and the fluid 
velocity. To define the geometry for the purposes of this paper, the 
complete set of geometrical data related to the heater oval shape, 
diameter, length and rod arrangement, as well as other minor 
geometrical details, may be approximately substituted by the local 
hydraulic radius. The density and the gas viscosity are included in the 
list of variables to take into account in calculating the inertia or form 
drag and the viscous effects, respectively. The pressure is included to 
consider compressibility effects that may be present at high veloci¬ 
ties. Thus, the following functional relationship can be applied: 


5p 

5x 


f(P,r he ,P,P,u) 


( 2 ) 


The integration of this function along the heater length leads to 
the instantaneous pressure drop at the heater A p. The process allows 
the gas flow to be characterised using an instantaneous, spatially- 
averaged velocity defined as u = m/(pA xxe ). Moreover, the list of 
influential parameters is completed with the gas temperatures at 
both the entry and exit sections, Tgo and T g L, as well as with the wall 
temperature, which usually is assumed to be uniform, to consider 
the spatial variations of gas properties. Alternatively, Tgo and T g i can 
be replaced by the gas temperature difference A T g and the mean gas 
temperature T g . Therefore, it can be written: 

Ap = f(p,R he ,p,p,u,L e ,AT g T g ,T w ) (3) 

The application of Buckingham’s theorem to the last expression 
(3) leads to the following functional relationship with six dimen¬ 
sionless groups, in which R/, e , p, u and T w have been used as refer¬ 
ence quantities: 

_Ap_ _ F ( t p Le A7 g 

1 _2 l P u Rhe 1 —2 Rhe ?w T w 

-pu z \ ne -pU z ne J 
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Table 2 

A comparison among the geometrical characteristic parameters for the prototype heater, the strictly similar derivative and the experimental heater. 



Vdxe (cm 3 ) 

Vde (cm 3 ) 

A wxe (cm 2 ) 

Rhe (mm) 

L e (mm) 

Axxe (cm 2 ) 

^dxe ( ) 

Mc/e ( ) 

a wxe ( ) 

(-) 

&xxe ( ) 

Rhel^e ( ) 

Ml 02C 

5.67 

16.77 

331.38 

0.171 

37.25 

1.52 

0.088 

0.261 

20.679 

0.0043 

0.0949 

0.0046 

Strict derivative 

23.22 

68.69 

848.17 

0.274 

59.60 

3.89 

0.088 

0.261 

20.679 

0.0043 

0.0949 

0.0046 

Relaxed derivative 

57.00 

76.85 

817.00 

0.698 

108.00 

5.28 

0.167 

0.225 

16.712 

0.0100 

0.1080 

0.0065 


The ideal gas model allows the pressure to be substituted by the 
product pRT re f, where T re / denotes the reference temperature at 
which the fluid properties are computed. Therefore, the pressure 
drop factor, sometimes called the Euler number, is a function of the 
Reynolds number, the Mach number, and the geometrical and 
temperature ratios: 


phenomena whose understanding depends on how the same 
dimensionless numbers are interrelated. 

For example, assuming isothermal conditions, Eqs. (4) and (5) 
reduce to the following well-known expressions that apply for 
the steady, fully-developed, laminar flow of incompressible, New¬ 
tonian fluids through smooth, straight ducts: 


jy - Fl N N Le ATg Tg 

— r \ l 'rei n mch~n ——if - 
Afte 1 w M 


(4) 


16 L t 
N^e R 


^P _ m 2 U_ 
P ma N re R he 


c = 16 
f N re 


The advantages of using the Stirling number to characterise the 
Stirling engine phenomena have been evidenced previously 
[21,26,27], but it is convenient to maintain the more common 
expression N re = N^alNsg in the last functional relationship to 
facilitate comparisons. 

A merely algebraic substitution allows the Eq. (4) to be modified 
to introduce the fractional pressure drop as follows: 


Ap 

P 


_ N z -Fl N N Le ^ Tg Tg 

— l ^ma 1 1 N re, iv ma , —, -=—, =- 
n he 1 w 1 w 


(5) 


An analogous procedure can be applied to analyse the instan¬ 
taneous convective heat transfer coefficient, which then requires 
that the list of independent variables of the Eq. (3) be completed 
with the addition of characteristic parameters of the gas thermal 
performance, namely, its specific heat capacities. Moreover, if the 
heat conduction in the heater rods is considered to influence the 
behaviour, the volumetric, specific heat capacity and the thermal 
diffusivity of the wall material must also be taken into account. The 
pressure is not explicitly included here because the product RT re f is 
determined by the gas temperatures and the specific heat capac¬ 
ities. Thus, the following expression can be written: 


h — f(Rh e , p, p, u, k, Cp, Cy, L e , AT g , Tg, T w , p w c w , a w ) (6) 

Applying Buckingham’s theorem to the last expression leads to 
the following functional relationship that involves ten dimension¬ 
less groups, in which Rhe, p, u and c p have been used as reference 
quantities: 


h p / (l Le k c v CpATg CpT g c p 7\v p w c w ol w \ 
pucp ~ \puRp e ’ R he ’ puc p R he ’ Cp ’ u 2 ’ U 2 ’ U 2 ’ PCp 'uRheJ 


Taking into account that c p /R = 7/(7 - 1) and manipulating the 
equation to obtain more meaningful dimensionless groups, this 
equation can then be written as follows: 


Na = F(N re ,^,N pr , 7 , 


A7g Tg 


N m l,J^,N TC R,Nf 0 


T w T, 


ref 


Thus, for T re f = T g , the previous equation is finally reduced to: 


N -Fl N N N r Le A7g Tg Pw c wAT g 0L W \ 

N st -F{N„,Npr,N ma , y.— 


(7) 


Eqs. (4), (5) and (7) are the basis of the experimental analysis 
described below. To correctly interpret the results, it must be noted 
that these equations can be applied to a variety of physical 


However, it can be noted that the introduction of Euler, Reynolds 
and Mach numbers in this problem is arbitrary because inertia, 
form drag and compressibility effects have not been considered. In 
other words, density and pressure are not influencing properties 
based on Hagen-Poiseuille’s law: 


6V 
5 x 



constant 


^=W L 

K he 


e 


The revision of the analysis after removing p and p from the list 
of variables shows that, in this case, only one dimensionless group 
is linked to the length-to-hydraulic radius ratio. That dimensionless 
group may be interpreted as a version of the Stirling number, which 
is the product of several of the abovementioned groups that, 
therefore, cannot vary independently: 


kpRhe = 
pu R he 



Neu Nre 



Ap 

P 


Nre-N 


-2 

ma 


Eqs. (4) and (5) can also be reduced to expressions that apply for 
cases where the density is needed for modelling. For example, 
Forchheimer’s law takes into account the pressure drop not only 
caused by viscosity but also by the acceleration and deceleration 
effects of the fluid as it travels through the tortuous flow path of the 
porous media: 


8p 

8x 


/l-+f}pu 2 

K 


{-dp/5x)/P 

pu 2 


p 

punfi 


+ 1 


The spatial integration of this law allows expressions to be 
obtained where the Reynolds number is independent either from 
the Euler number, the product between the fractional pressure drop 
and the inverse square of the Mach number, or from the friction 
factor, that is to say: 


Neu 


U + b ) Rhe 

- N 2 f— + b\ 

2 ma \N re + ) 



L e _ 

Rhe 




r =-°- 

f N rP 


+ b 


In short, both Hagen-Poiseuille’s and Forchheimer’s laws are 
models of incompressible flows in which Ap/p is an explicit function 
of the squared Mach number and C/ does not depend on N ma . 

However, Eqs. (4) and (5) also agree with the experimental 
measurements carried out by Su [28] for steady, compressible flow 
through stacked screens, as well as with the analyses for Stirling 
regenerators introduced by Organ [ 16,17,27], This author recommends 
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that the experimental correlations for regenerators be presented using 
N ma as a parameter because it is influential at values as low as 
N ma ~ 0.1. Although the compressibility effects are not probable at 
moderate velocities, Organ’s recommendation will be followed in the 
present paper to show how those effects can be evaluated through the 
different functional influences of N ma , as seen for the aforementioned 
models of incompressible flows. Therefore, the independence of both 
Cf and N st on N ma , as well as the dependence of Ap/p on N^a. are 
expected in the non-tubular heater tests. 

4.1 Experimental test rig 

The application of friction factor and Stanton or Nusselt number 
steady flow correlations for Stirling engine heat exchangers is open 
to discussion because the Stirling engine operation implies not- 
fully-developed, bidirectional flow with a variable mass rate. 

For example, Gedeon and Wood [29] pointed out that there is no 
essential difference between stationary and oscillating flow under 
the conditions used to test a few regenerators. Conversely, Isshiki 
et al. [30] concluded that friction factors for oscillatory flow are 
larger than those in steady flow. In the same way, Ibrahim et al. [31 ] 
studied the pressure drop through complex geometries under 
oscillatory flow conditions and found that the experimental results 
for oscillatory flow were 1.2—1.3 times the predicted values for the 
steady operation. 

Regardless, the consensus is that steady flow tests may be 
meaningful to characterise Stirling engine heat exchangers. 
Therefore, an experimental test rig was constructed (see scheme in 
Fig. 4) to obtain the correlations for the experimental heater based 
on the steady flow adaptations of Eq. (4) or (5) and (7). 

The compressed air used in the tests was pre-processed at the 
supply network for moisture prevention and filtering. An additional 
filter was also arranged at the entry of the mass flow meter, which 
allows the mass flow to be controlled in the range from 0 to 19 kg/h, 
with variations of 0.02 kg/h. 

The absolute air pressure was measured at the heater entry 
using a transducer controlled in the range of 1—6 bar through two 
valves, the first one placed at the connection to the compressed air 
network and the other one, at the heater exit. The pressure drop 
across the heater was measured using a water column differential 
manometer that was also controlled by those valves in the range 
of 0—1200 Pa. The full-scale accuracy is ±0.25% for the pressure 
transducer, and ±0.5% for the pressure drop gauge. 

The air temperature was calculated at both the entry and exit 
heater sections and the average value of four measurements was 
obtained by type I< thermocouples (200 pm wire diameter). The 
wall temperature was measured as close as possible to the bases of 


the heater rods and to the heater centre using other similar ther¬ 
mocouples inserted in a hole. An additional thermocouple provided 
the feedback signal for the electrical resistance controller used for 
heating. The thermocouples used were manufactured under the IEC 
584-2 (1982) Standard, which implies an accuracy of ±1.6 °C for the 
temperature range measured. 

4.2. Experimental results 

The experiments were carried out at heater wall temperatures 
of 100 °C, 200 °C, 300 °C and 400 °C. The gas pressure and mass 
flow rate were regulated to obtain series of constant N ma from 
0.0010 to 0.0180, which implies values of N re between 60 and 1400. 
The temperature ratios oscillated between 0.045 and 0.325 for A T g l 
T w and between 0.725 and 0.950 for T g /T w . Due to variations in the 
gas properties, N pr varied in the range from 0.745 to 0.765. 

Fig. 5 shows the fractional pressure drop values obtained for 183 
experimental points measured in the abovementioned ranges. The 
following correlation is a particular case of the Eq. (5) for the 
experimental heater value of Rhelle = 0.0065 and adjusts to the 
measurements with an R-squared value of 0.9985 and an RMSE of 
7.81%: 

An /T \ ~°- 9 41 /AT \ -0.033 

f = 174.058NL 079 N re 0 ' 238 g) (8) 

Fig. 5 must be interpreted by first noting that data with the same 
N m a symbol, in general, have different values of temperature ratios. 
Moreover, the lines of constant N ma were calculated from Eq. (8) for 
the mean values of the temperature ratios, that is to say, for 
Tg/T w = 0.86 and AT g /T w = 0.200, to provide some graphical idea 
of the dispersion caused by these ratios. 

Conversely, the experimental N ma exponent obtained, which 
was close to 2, and the inspection of Eq. (5), seem to confirm that 
the compressibility effects are negligible. Additionally, the near¬ 
zero exponent of A T g /T w reflects its limited influence on the frac¬ 
tional pressure drop. Thus, taking these pieces of evidence together 
suggests using the following friction factor correlation, which is 
simpler than Eq. (8) and has sufficient accuracy in practise because 
it adjusts to the experimental data with an R-squared value of 
0.9990 and an RMSE of 9.02%: 

S -0.991 

Cf = 0.904N re ° 145 LM (9) 

Fig. 6 shows the experimental data and prediction lines of the 
constant T g /T w deduced from this equation. 
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Fig. 4. The schematic of the experimental test rig. 
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Fig. 5. The fractional pressure drop data obtained from experimental tests. 


Figs. 7 and 8 show the Stanton number values obtained for the 
183 experimental heat transfer points measured in the above- 
mentioned ranges. The following correlation is a particular case of 
Eq. (7) for the experimental heater value of RhelL e = 0.0065 and 
adjusts to the measurements with an R-squared value of 0.9958 and 
an RMSE of 7.47%: 

/at \ 0.856 \ 4.940 

Nst = 0.0176< 091 JV m r 27 Ni r 870 ^J g) (10) 

The steady flow conditions justify the assumption that the 
influences of the transient heat conduction variables, that is to say 
N T cr and N/ 0 , were not noticeable. 

Fig. 7 shows the lines of constant Mach number deduced from 
the correlation for the constant values T g /T w = 0.86, A T g l 
T m = 0.200 and N pr = 0.75. Fig. 8 shows the lines of constant A T g /T w 
deduced from the correlation for the constant values of N pr = 0.75, 
N M a = 0.0075 and T g /T w = 0.86. In this situation, N st increases 


slightly with N re , while N ma has little influence, and the dispersion 
of the results is caused mainly by the temperature ratios AT s /T w and 
T g /T w . 

Therefore, as seen in the fractional pressure drop, the 
compressibility effects associated with the Mach number influence 
seem negligible, which suggests using the following simpler 
correlation that adjusts to the experimental data with an R-squared 
value of 0.9955 and an RMSE of 7.57%: 

/at \ 0.856 x 4.975 

Nst = 0.023JV l ° e ° 66 JVp r 750 (jAj (11) 

Additionally, unlike in tubular geometries, in this situation, the 
Reynolds analogy is not fulfilled. 

Conversely, the comparisons with classical correlations show 
the considerable influence of the temperature ratios in Eqs. (9) and 
(11). Pending a qualitative explanation that might be derived from 
further research, these tests have conclusively shown the need to 



REYNOLDS NUMBER 


Fig. 6. The comparison between the experimental friction factor and the prediction lines for a constant T g /T w deduced from Eq. (9). 
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Fig. 7. The comparison between the experimental Stanton number and the lines of constant Mach number deduced from Eq. (10) for constant values of N pr = 0.75, T g /T w = 0.86 
and \T t ,/T w -0.200. 


include the temperature ratios as correlation parameters to obtain 
accurate results, and the decrease in accuracy if the exponents of 
those ratios are diminished. Fig. 9 shows that calculations based on 
the Eq. (11) fit most experimental data within to ±10% of the 
dashed lines, while, for example, the following equation adjusts to 
the measurements with an R-squared value of 0.9818 and an RMSE 
of 17.5%, with most experimental data within ±25%: 

Nst = 0.00086< 391 Ni r 5 ^j (jj) 05 (12) 

5. Discussion 

A complete analysis of the predictions of Eqs. (9) and (11) from 
the viewpoint of the engine operation is outside of the scope of this 
paper. However, Table 3 provides some approximate idea by 


comparing these predictions with the friction factor and the Stan¬ 
ton number values of benchmark prototypes. 

Turbulent flow is usually present in Stirling engine heaters [32], 
and the Reynolds number reaches its maximum values of around 
5100 in the M102C engine heater at the operating point assumed 
for scaling (Table 1 ), and around 9400 in the heater of the relaxed 
derivative, due to its out-to-scale, higher swept volume. The rest of 
the Nre.max values in Table 3 have been estimated from the works of 
Seume and Simon [32] and Fano [33]. The correlation assumed to 
calculate the friction factor for slotted or tubular heaters has been 
recently used by Organ [34] to evaluate the available work lost due 
to pumping losses in Stirling engine heat exchangers, namely 
C f = 0.462 Nre 02 . The classical correlation by Dittus and Boelter was 
used for the Stanton number, namely N s t = 0.023 Nr e 0 2 NpP £67 . 

The comparisons show that the friction factor values are higher 
for the experimental heater than for the tubular heaters in other 
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Fig. 8. The comparison between the experimental Stanton number and the lines of constant AT g /T w deduced from Eq. (10) for constant values of Nma = 0.0075, N pr = 0.75 and 
Tg/T w — 0.86. 
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Correlation Stanton number 


Fig. 9. The comparison between the experimental Stanton number and the correlation estimations based on Eq. (11). 


Table 3 

A comparison between the friction factors and Stanton numbers for heaters of 
benchmark engines. 


Prototype 

Nre, max 

Friction factor, Cf 

Stanton number, N st 

4-215 

-200,000 

0.04 

0.002 

P-40 

-100,000 

0.05 

0.003 

GPU-3 

-10,000 

0.07 

0.004 

M102C 

-5100 

0.08 

0.005 

Relaxed derivative 

-9400 

0.28 

0.003 


engines, which probably cancels out the aforementioned advantage 
of having a favourable RhelL e ratio. Conversely, the Stanton number 
estimations provided by Eq. (11) seem more acceptable, and it 
seems interesting to note that correlations with lower exponents of 
temperature ratios, like Eq. (12), would lead to higher N s t values. 

In short, from the standpoint of both pressure loss and heat 
transfer, it seems that a heater with the type of geometry proposed 
cannot overcome the characteristics of a tubular model. However, it 
cannot be excluded that an optimised non-tubular heater could be 
interesting for applications where high Reynolds numbers are 
reached in the heater and there are practical problems associated 
with the use of bundles of tubes, given that the available work lost 
due to heat transfer irreversibilities in tubular heat exchangers 
usually overwhelms the pumping loss by an order of magnitude 
[34], 

Perhaps the most advisable alternative to optimise the experi¬ 
mental heater geometry is to validate a 3D-CFD model through 
experimental results, like those provided by this paper. The CFD 
technique is finding increasing application to the modelling of both 
Stirling engine components [35,36] as the global engine perfor¬ 
mance [37]. Such a model could be the basis for analysing the range 
of performance for which the non-tubular geometry is appealing, 
taking into account the reciprocating character of the gas flow in 
the actual engine operation. 

6. Conclusions 

A micro co-generation solar unit has been designed by scaling 
the Philips M102C beta prototype, assuming a case of relaxed 
similarity. 


The dead volume obtained in the non-tubular approach 
considered for the derivative heater is more than twice the value of 
the dead volume in the strictly similar heater. This disadvantage has 
been compensated for by the reduction of dead volume in the 
adjacent spaces, which resulted in a computed value £o ~ 0.302 for 
the relaxed derivative, which is slightly higher than the value for 
the prototype, despite the different engine configuration and drive 
mechanisms. 

The experimental heater ratio RhelLe is 41% higher than the 
strictly similar value, but the dimensionless characteristic param¬ 
eter of heat exchange a wxe is 19% lower for the relaxed approach, 
despite the fact that the wetted surface is almost identical in both 
cases. 

The variables needed to characterise the friction and heat 
transfer performance of the experimental heater have been iden¬ 
tified, appropriate functional relationships have been derived by 
means of dimensional analysis, and experimental data have been 
obtained in steady flow conditions. 

The N ma exponent, which is close to 2 in the fractional pressure 
drop correlation, seems to confirm that the compressibility effects 
are negligible. Additionally, the near-zero exponent of AT g /T w 
reflects its limited influence on the fractional pressure drop. 

The following friction factor correlation adjusts to the experi¬ 
mental data with an R-squared value of 0.9990 and an RMSE of 
9.02%: 

/T s -0.991 

Cf = 0.904JV re 0145 (=«J (13) 

The compressibility effects associated with the Mach number’s 
influence seem also to be negligible in the Stanton number 
measurements, which suggests using the following correlation that 
adjusts to the experimental data with an R-squared value of 0.9955 
and an RMSE of 7.57%: 

/AT \ 0 856 /T \ 4.975 

Na = 0.023< 066 iVi r 750 ^J (-A (14) 

The Reynolds analogy is not fulfilled by the experimental heater, 
and the friction factor values predicted by the Eq. (13) are higher 
than those estimated for the tubular heaters of benchmark engines, 
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which probably cancels out the aforementioned advantage of [18 
a favourable Rh e /L e ratio. Conversely, the Stanton number estima¬ 
tions seem more favourable. [ ig 

An optimised non-tubular heater may still be interesting for 
applications where high Reynolds numbers are reached in the 
heater and there are practical problems associated with the use of 
bundles of tubes. An advisable alternative to optimise the experi¬ 
mental heater geometry would be the construction of a 3D-CFD I 21 
model, validated with abovementioned experimental data. 

[22 
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